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The gap equation for Dirac quasiparticles in monolayer graphene in constant magnetic and pseu- 
domagnetic fields, where the latter is due to strain, is studied in a low-energy effective model with 
contact interactions. Analyzing solutions of the gap equation, the phase diagram of the system 
in the plane of pseudomagnetic and parallel magnetic fields is obtained in the approximation of 
the lowest Landau level. The three quantum Hall states, ferromagnetic, antiferromagnetic, and 
canted antiferromagnetic, are realized in different regions of the phase diagram. It is found that the 
structure of the phase diagram is sensitive to signs and values of certain four-fermion interaction 
couplings which break the approximate spin-value 51/(4) symmetry of the model. 


I. INTRODUCTION 

Among many remarkable properties of graphene its response to strain can be singled out as one of efficient means 
to change and control the characteristics of the electronic states in graphene. Since graphene is only one atom thick, it 
is easily subjected to mechanical deformations. Various proposals to engineer strain in graphene were discussed in the 
literaturei^— . It is known that strain induces effective gauge fields^i^ with the corresponding effective ’’magnetic” fields 
of opposite sign in valleys K and K' that means that elastic deformations, unlike real magnetic field, preserve time 
reversal symmetry^ii. Since time reversal symmetry is unbroken, strain induced fields are known as pseudomagnetic 
fields in the literature (for a review of gauge fields in strained graphene, see Refs. Hi). It was proposed in Refs. [131111 
that a designed strain may induce uniform pseudomagnetic field, which can easily reach values exceeding 10 T. 

The observation of anomalous integer quantum Hall (QH) effect with the hlling factors v = ±4(|n| -|- 1/2) ( n is the 
Landau level index) in graphene in a magnetic fieldi^, in accordance with theoretical studies in Refs. |l3lll4l . was a 
milestone in graphene research as it became a direct experimental proof of the existence of gapless Dirac quasiparticles 
in graphene. The four-fold degeneracy of the Landau levels in graphene is due to the 517(4) symmetry connected 
with valley and spin. Later the plateaus i/ = 0,±1,±4 in the QH effect in graphene were observe d^^d^ in a strong 
magnetic held B > 20T. These plateaus are connected with the magnetic held induced splitting of the n = 0 and 
n = 1 Landau levels and the degeneracy of the lowest Landau level (LLL) is thus completely resolved. 

The Landau levels related to pseudomagnetic heldsiSiii were observed in spectroscopic measurements It was 
pointed out^^ that pseudomagnetic helds due to strain can interfere in many ways with real magnetic helds. For 
example, the interplay of pseudomagnetic and magnetic helds in the quantum Hall regime causes backscattering in 
the chiral edge channels that can destroy the quantized conductance plateaus and gives rise to unconventional QH 
ehect in strained graphenei^ with oscillating Hall conductivity. 

The gap generation in graphene in the presence of a pseudomagnetic held was studied in Ref. . Interestingly, 
it was found that unlike magnetic held which catalyses the generation of the time reversal invariant Dirac mass, 
pseudomagnetic held catalyses the generation of time-reversal symmetry breaking Haldane mass. Various competing 
ground states in monolayer graphene in pseudomagnetic helds were recently studied in Ref. [2lj| . Finally, we would 
like to add that very strong 50 — 60 T pseudomagnetic helds may be realized in molecular graphene^^. 

The interplay between different possible ground states in strained graphene in a magnetic held represents an 
important unsolved problem at the moment. In the present paper, we study a gap generation for quasiparticles in 
monolayer graphene in the presence of both constant magnetic and pseudomagnetic helds and using the model with 
local four-fermion interactions considered in Refs. f^l^ . Local four-fermion terms in the Hamiltonian are remnants 
of the interactions on the atomic scale, and in spite being much smaller than Coulomb interaction, they play an 
important role in deciding how the 5?7(4) symmetry is broken in monolayer graphene as well as bilayer graphene. 
This especially concerns the nature of the QH state with half-hlled zero-energy Landau level. We obtain the phase 
diagram for competing quantum Hall states in the LLL approximation when the chemical potential is tuned to the 
charge neutrality point, i.e., the state with the zero filling factor. 

The paper is organized as follows. We begin by presenting in SeclUthe model describing low-energy quasiparticles 
excitations in strained monolayer graphene in an external magnetic field and in the presence of local four-fermion 
interactions. The derivation of the gap equation is given in Sec IHII and its solutions are presented in Sec lIVI The phase 
diagram of the system is derived and discussed in SeclV] The main results are summarized in Sec IVIl Appendices at 
the end of the paper contain technical details and derivations used to supplement the presentation in the main text. 
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II. HAMILTONIAN OF THE MODEL 

The low-energy quasiparticles excitations in graphene can be described in terms of a four-component Dirac spinor 
'^a — {'^KAaT'^KBa,4’K'Ba,4’K'Aa) which combines the Bloch states with spin index a = 1,2 on the two different 
sublattices (A, B) and with momenta near the two nonequivalent valley points {K, K') of the Brillouin zone. The free 
quasiparticle Hamiltonian has a relativistic-like form with the Fermi velocity vp = 10® m/s playing the role of the 
speed of light 


Hq = J + 7 ^%)^ + , ( 1 ) 

where ^ = vj/t-yO is the Dirac conjugated spinor and r = {x,y). The matrices 'j" with f = 0,1,2 are 4x4 matrices 

which satisfy the anticommutation relations of the Dirac algebra { 7 '', 7 “^} = , where = diag(l, —1, —1) and 

/i, F = 0,1, 2. These matrices belong to a reducible representation of the Dirac algebra = fz® {Tz,iTy, —iTx), where 
the Pauli matrices fi and Ti with i = x,y,z act in the subspaces of the valley {K, K') and sublattices (A, B) indices, 
respectively. 

The canonical momentum tt = —i/iV -I- cKjc -V y^y^eAs/c includes the vector potential in the Landau gauge 
A = (0, B±x) corresponding to the component Bj^ of an external magnetic field B orthogonal to the plane of 
graphene, and A 5 is the vector potential describing the strain induced gauge field o^®i^°'^®i^® . In the representation of 
the Dirac matrices that we use, 7 ^ and 7 ® = matrices equal y® = ify ® I and = fx ® I, where I is the 

2x2 unit matrix. It is easy to check that the product of these matrices y®y® = 0 / is a matrix diagonal in the 

subspace of valleys that ensures that the term in the canonical momentum with the vector potential A 5 takes opposite 
signs in K and K' valleys. In what follows, we will consider only the case of constant magnetic B and pseudomagnetic 
B 5 fields. The pseudomagnetic field B 5 points always in the direction perpendicular to the plane of graphene and, 
therefore, is described by the vector potential A 5 = (O^B^x), where B 5 = jBsj. 

The last term in the free Hamiltonian ([T]) is the Zeeman interaction ez = /tbB with g,B = e?i/ (2mc) being the Bohr 
magneton and az is Pauli spin matrix whose eigenstates describe spin states directed along or against the magnetic 
field B. Here B = + B^| is the strength of the magnetic field and B|| is its component parallel to the plane of 

graphene. We note that the standard Zeeman interaction ^^Bcr can be reduced to this form using a rotation in spin 
space. 

The Coulomb interaction between electrons is described by the following Hamiltonian: 

He = — J (Prd'^r’^{r)y^'9(r)Uc{Y — 

where Uc{x) is the Coulomb potential. In order to simplify the analysis, we follow the approach of Ref. and replace 
the Coulomb interaction C/c(r) by the contact interaction Gmt^^(r). The Hamiltonian Hq + He in the absence of the 
Zeeman term possesses a global 51/(4) symmetry connected with valley and spin degrees of freedom. 

Although the Coulomb interaction is the strongest interaction between electrons in graphene, local four-fermion 
interaction o^^i^® play a crucial role too. Although these interactions are much smaller than the Coulomb one, they 
break, in general, the 5C/(4) symmetry and crucially affect the selection of the ground state of the system. A set 
of local valley and sublattice asymmetric four-fermion interactions was introduced in Ref. [ 2 ^ . The 1 / = 0 quantum 
Hall state was studied and it was shown that the phase diagram, obtained in the presence of generic valley and 
sublattice anisotropy and the Zeeman interaction, consists of four phases: ferromagnetic, canted antiferromagnetic 
(CAF), charge density wave, and Kekule distortion. The Hamiltonian of generic local four-fermion interactions reads 

dicontact ^ J (fr^gjk[i{r)y°Tjk'^{r)]'^, Tjk=fj®Tk, ( 2 ) 

where j, k = x,y, z. We do not include in Hcontact the term with goQ as it corresponds to the local Coulomb interaction, 
which has already been taken into account by Gmt- In addition, we dot not include in our model the terms with gofc 
and gjo, which vanish in the first order in the Coulomb interactions and arise only in the second order due to virtual 
transitions to other bands^. The coupling constants gjk are not all independent. As shown in Ref. [ 2 ^, symmetry 
and other considerations lead to the following equalities for nonzero constants: 

= 9xx = 9xy = 9yx = 9yyi 9-Lz = 9xz = 9yz^ 9zX = 9zx = 9zy ( 3 ) 

Thus, totally we have four interaction coupling constants, Gint, 9±±, 9±z, 9z±j in the considered model. Finally, let 
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us present Tjk in terms of the 7 -matrices 

%x = -*7^7^, %y = *7^7\ Txz = - 7 ^ 7 °, 

Tyx = - 7 V, Tyy=-f^j\ Tyz=il^l°, 

Tzx = Tzy = -i7\ Tzz = 7°- (4) 

All these matrices are normalized as = 1. This presentation is useful for the derivation of the gap equation in the 
next section. 


III. GAP EQUATION 

We will solve the gap equation in the Hartree-Fock (mean-field) approximation^— which is conventional and 
appropriate in this case. In the subsection IIII Al we will derive the gap equation in the case where only real magnetic 
field is present. In the next subsection, we will generalize the gap equation to the case where both magnetic and 
pseudomagnetic fields are present. 


A. Magnetic field 


At zero temperature and in the clean limit (no impurities), the Schwinger-Dyson equation for the quasiparticle 
propagator G{u,u') = /i“^(0|T4'(M)^(u')|0) in graphene in the mean-field approximation takes the form 


iG ^{u,u') = iS — hGintj'^G{u,u)j^S{u — u') + hGint^^ti [j'^G{u,u)] S{u — u') 

- h'^gjk{j°TjkG{u,u)-f°Tjk - 'y°Tjktr [y°TjkG{u,u)]}S{u - u'), 


( 5 ) 


j,k 


where u = (t, r). In this subsection, we will derive the gap equation in graphene in a magnetic field. The generalization 
to the case of both magnetic and pseudomagnetic fields is rather straightforward and will be considered in the next 
subsection. 

The inverse free propagator in the case under consideration is given by 

u') = [{ihdt - ez<7z)l° - vf{t^ ■ 7 )]Siu - u'). ( 6 ) 


For the full quasiparticle propagator, we will use an ansatz which is a generalization of the ansatz used in the previous 
work by two of us 21 

iG~^{u,u') = + Ai 7 ° 7 ^ 7 ^ — u_F( 7 r • 7 ) — A -|- S{u — u'), (7) 


where matrices /r,/i. A, A are defined as /r = = A^a,j,A = A^u^, and index v runs the values 

V = 0,x,z with ax and az being Pauli spin matrices and ao the unit 2x2 matrix [the absence of quantities with 
ay matrix is consistent with subsequent analysis of a gap equation]. In what follows, we consider twelve dynamically 
generated parameters /Tjy, fl^, A^, and Aiy as constant that is consistent with our mean-field analysis of the present 
model with contact interactions. 

The parameters fij and flj with j = z are generalized chemical potentials connected with the QH 
ferromagnetis m^^i^^~ — . On the other hand, Aj and Aj are related to the magnetic catalysis scenario^^— and 
are Haldane and Dirac masses, respectively, and correspond to excitonic condensates (for a brief review of the QH 
ferromagnetism and magnetic catalysis scenario, see Refs. [39l - l4lj |). Actually, it was shown in Ref.jl^ that the QH 
ferromagnetism and magnetic catalysis scenario order parameters necessarily coexist. The physics underlying their 
coexistence is specific for the systems with relativistic-like quasiparticle spectrum that makes the quantum Hall dynam¬ 
ics of the SU (4) breakdown in graphene to be quite different from that in conventional systems with non-relativistic 
quasiparticle spectrum. 

According to Eq.(I7]), the full propagator G{u,u') can be written in the form 


G(m, u') = i{u\ {ihdt + — vf{tv ■ 7 ) -|- i/ly^y^ -|- iAy^y^y^ — A jit') 


-1 


( 8 ) 


where the states |u) are eigenstates of the time-position operator u: u\u) = u\u), {u\u') = 5{u — u'). In Appendix (Al 
we derive an explicit expression for the propagator G{u,u') in the form of a sum over Landau levels. 







4 


The symmetry-breaking generalized chemical potentials and gaps Aj^, are related to the corresponding 

order parameters through the following relationship: 

=-htr[OM{u,u)], (9) 

where 8x8 matrices Oi, = 7 °crjy, ctj,, respectively. Compared to our previous analysis j^^d^ we 
included the spin matrix in order to be able to describe the canted antiferromagnetic state. 

Since the right-hand side of the gap equation ([S]) contains the full propagator at the coincidence limit u' = u, we 
should calculate G{u,u')\u=u' = G{u,u), where G is the translation invariant part of the full propagator defined in 
the mixed frequency-momentum representation in Ea. (IA15|) . By making use of Eqs. (IA19|) and (IA20I) . we find that 


G{u, u) 



dujd?k 

(27r)3fi, 


G(a;,k) 


i \ ^ f duj [P- -I- P+9{n - 1)] 

2TrP J 2Trh A4 — ne\ 


( 10 ) 


where I = ^/hcJ\eB±\ is the magnetic length, es = ^j2hvp\eBj_\/c ~ 424A/|i?_L|[T]K is the Landau energy scale, P± 
are projectors given by Ea. IjAlip . and W, M are matrices expressed through /i, jl, A, A and defined in Eas. (IA2l) and 
(IA3I) . We note that the filling factor iz = 2'Kpp is related to the carrier imbalance p = Ue — Uh, where Ue and Uh are 
the densities of electrons and holes, respectively, and p is determined through the Green’s function as 

p = (0|'I't(u)'l>(M)|0) = —?itr[y°G(M,'u)]. (11) 

Since W and A4 contain only y*^ and y^y^ Dirac matrices, it is convenient to work with eigenvectors of these 
matrices. The equality (y^y^)^ = —1 implies that the eigenvectors |si 2 ) of the matrix y^y^ are purely imaginary 

y^y^|si 2 ) = *Sl 2 |si 2 ), Sl 2 = ±1. (12) 

Similarly, since (y°)^ = 1, the eigenvectors of the matrix y° are real and given by 

y” I So) = So I So), so = ±l. (13) 

Furthermore, since y*^ and y^y^ commute, we can consider states |sosi 2 ) which are simultaneously eigenvectors of y*^ 
and y^y^ with eigenvalues Sq and *Si 2 , respectively. The vectors |soSi 2 ) form a complete basis, and since W and A4 
contain only y° and y^y^ matrices, the propagator G(u,u) is diagonal in the basis of |si 2 So) vectors and is given by 

G(m,m;so,Si 2 ) = J + m) + d][co + 

X (l + Si 2 S±-b [1 - Si 2 S_L] 6 >(n-1/2)^, j = x,y,z. (14) 

Here s_l = sgn(ei?j^), the matrices m and d are defined in Appendix [Al and the coefficients co,Cj are 

Co = u!^ + 2ujmo — ne% + m'l — dl, Cy = 2isQ{dxmz — d^nix), (15) 

= 2{ujm^ + moTTix - dodx), = 2{u}m^ + mom^ - dod^), (16) 

where = pv — soSi 2 Pv, di, = A^, -b S 0 S 12 A 1 / (z/ = 0, x, z) and summation over dummy index v is meant. For strong 
magnetic fields, we write 


G{u,u) = Glll{u,u) + GhLL{u,u), (17) 

where we separated the contributions of the zero Landau level, Glll{u,u) with n = 0, and higher Landau levels, 
GhLL{u,u) with n > 1, in Ea. dMl) . 

Let us calculate first the lowest Landau level propagator Glll{u,u). In order to integrate over uj in Eg. lfni) . we 
rewrite the integrand by using the relation 


Co -b Cjaj = [so(w -b m) — d] [so(a; -b to) -b d] 


(18) 


valid for n = 0 and assume as usual that uj is replaced by w -b iesgnw and e —?► 0+. Hence we obtain that Ea. dTH) 
implies the following propagator at the limit of coinciding points in the LLL approximation: 


Glll(u,u; So, 312) = 


i l-bSl2S_L 
^ 2 


— CX) 

So 1 -b Si 2 S L 


r du, ^ 

{rrii - sodi)a, 


0{E- \pl\) + sgn{pL) d{\pL \ - E) 


AnhP 2 


E 


i = x,z (19) 
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where the factor (1 + si2S_l)/2 with s_l = sgn(ei3j_) reflects the presence of the spin projector P_ = (1 — fs_L7^7^)/2 
in the LLL contribution, energy 


E = - s^d^Y + (m^ - So4)2, 


and we introduced the notation 


— fno — sodo — fiQ — soSi 2 /io ~ sqAo — S12A0 


( 20 ) 


( 21 ) 


for an ’’effective chemical potential” in the lowest Landau level. 

By integrating over w, it is not difficult to check that the higher Landau level contribution GuLhiu^u) diverges as 
. Indeed, making the change of the variable lo —>■ and taking into account that all the dynamically 

generated parameters are much less than the scale es, we find that the leading contribution at large ne\ is given by 


GhLE^U, u; So, S12) 


1 


E 


1 djyO'i/ 
y/n tB 


( 22 ) 


This means that the right-hand side of the gap equation ([5]) diverges too. This result is the well-known artefact of 
using a model with local four-fermion interactions. For a long-range interaction like, for example, the Coulomb one 
considered in Ref.jd^, such a divergence is absent because the gap equation contains the quasiparticle propagator 
at different points u and u'. To proceed further, we regularize the divergence in the model under consideration 
introducing a cutoff nmax in the sum over Landau levels (a slightly different approach was used in Ref.fl^), which is 
connected with the ultraviolet (UV) cut-off in energy A (band width) according to the relation rimax = A^/e^. By 
using the regularization described above and retaining only the leading contribution, we find that the higher Landau 
levels contribution to the propagator at the limit of coinciding points is given by 


GhLLiu, u; So, S 12 ) 




(23) 


By combining Eas. (l23)) and (IT^ . the gap equation (O takes the following final form: 


SoTTlj^Ci/ di,(Tu ~\- SoSj_6z^z 


(Gint + ffzz) G(so, S 12 ) — ‘ig±± G(so, —S 12 ) + 2gz± G(—so, —S 12 ) 

2g±z G(—so, S 12 ) + gzz Y. (24) 

Sq,sJ 2 


where 


G(so,Si 2 ) = GlLl{u,U-,So,Si 2 ) +GhLL{u,U-,So,Si 2 ) 


-So 

AttHE 


{rrii - Sodi)(Ji 
E 


- IulD + sgn (^i) 0{\gL\ - E) 


1 + Si2S L 


A 


A-kH^v^p 


du Ou 


(25) 


and trace in the last term in Eg. (1241) is taken over the Pauli spin matrices [note that the quantities E and g,L depend 
on So, S 12 according to Eas. (l20)) . (|2T|) and rriu, d^ depend on them too]. In deriving Ea. (l24)l . we omitted the third term 
on the right-hand side of Eq.(IS]) which defines the Hartree contribution due to the charge density of carriers. The 
point is that there are other contributions due to the charges of ions in graphene and the charges in the substrate and 
gates. In view of the overall neutrality of the system, all these contributions should cancel exactly (the Gauss law). 

Einally, we would like to note that it is advantageous in deriving Ea. dMl) to use Tjk matrices given by Eq.(|3|) and 
utilize the Dirac algebra in order to calculate the contribution to the gap equation due to the last term in Eq.(IS|). 
[The matrices Tjk have simple commutation relations with the Green’s function G{u, u) which contains only 7 ° and 
7 ^ 7 ^ matrices (see Ea. (fTUl) L] 


B. Magnetic and pseudomagnetic fields 

In this subsection, we will derive the gap equation for quasiparticles in graphene in the case where both constant 
magnetic and pseudomagnetic fields are present. Then the canonical momentum has the form tt = —ifiV -I- eA/c -I- 
7 ^ 7 ®eA 5 /c, where A = (0,il^a;) and A 5 = {0,Box). Repeating the same computations as in the previous subsection, 
one can show that the only difference between the former and present cases is that the magnetic field B± is now 
replaced by the effective field B± + (i 7 ° 7 ^ 7 ^)i ?5 or B± — S 0 S 12 B 5 in the eigenstate basis of the matrices 7 ° and 7 ^ 7 ^. 
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Consequently, the pseudoinagnetic field has opposite signs in the K and K' valleys. Therefore, in order to take into 
account pseudomagnetic field, we should simply make the replacement B± —>■ Bj^ — S 0 S 12 -S 5 in the corresponding 
equations of Subsec ITlI Al except the Zeeman energy, where ez = ^b ; which includes the component of 

the magnetic field parallel to the plane of graphene. We find it convenient in the analysis below to use the notation 
&II = S||/Bj_ and 65 = B^/B^. 

Taking into account all Landau levels contributions, the gap equation for quasiparticles in graphene in the presence 
of constant magnetic and pseudomagnetic fields is given by 

som^ai, - + sos^ezcr^: = - {Gint + gzz) G^^\so,si2) - 4 g_L_L (sq, -S12) 

+ 2g,iG(5)(-so,-si2) + 25^,G(5)(-so,si2)+5.. E tr G^^) (s', ^), (26) 


where 


G(®nso,si2) = 


-so|ei?j_ - 5o5i2ei?5| [ {rrii - sodi)a^ 
4:TTh?C E 

1 + si2Sgn(ei?j_ - soSi2ei35) 


e{E- I^lI)+ sgn(^L) - E) 


A 


A-kH^v^p 


di/O u . 


(27) 


The first term on the right-hand side of Eq. (071) is the LLL contribution and the last one describes the contribution 
due to higher Landau levels. 

We have found solutions of the gap equation taking into account the contributions due to all Landau levels. It 
turned out that the higher Landau levels contribution in the weak coupling regime does not qualitatively change the 
results obtained in strong magnetic fields when the LLL approximation is valid. On the other hand, the higher Landau 
levels contribution essentially enlarges formulas and makes them very complicated. Therefore, in what follows, we will 
solve the gap equation and present our analysis in the LLL approximation omitting the contribution due to higher 
Landau levels. 


IV. SOLUTIONS OF GAP EQUATION IN THE LLL APPROXIMATION 

In this section, we consider solutions of the gap equation (I26p retaining only the LLL contribution. Propagator 
(E3 in the LLL approximation contains different projectors depending on which field i?j_ or B^ is stronger. We will 
consider both possibilities separately. 


A. \Bp\> IBsl 


Let us find solutions of the gap equation (1261) in the LLL approximation in the case where magnetic field is stronger 
than pseudomagnetic field. We will use the notation of Refs. [^l^ 


G,> 


Uo = 


2 ttP^ 


Uz = 


2 ttP^ 


u± = 


gj-z 

2 t:P' 


(28) 


The gap equation in the LLL approximation is obtained from Eg. (1261) by replacing the full fermion Green’s function 
G^®^ with G^ppp and multiplying both sides by the LLL projector (I-l-si 2 S ±)/2 (for specificity, in what follows, we take 
sx = +)■ It is easy to see that in the LLL approximation twelve parameters fj,^, and enter in combinations 

— A^) — {jli, — A^). Therefore, we have only six independent variables /i^ — Ai, and jlv — A^. Without loss of 
generality we can put fj,^ = fl^ = 0 so that rrii, = 0 and we are left only with parameters d^. Then the gap equation 
takes the following form: 


- d^{so)(7,y + soezCTz = - {^-sob5)-iuo 


+ Uz) 


di{so)<Ti 

E{so) 


0{E{so) - Mo(so)l) + sgn(do(so))6'(Mo(so)| - ^^(so)) 


-I- (1 -I- Sob5)uj_ 


di{ ^o)u'i 


g(E(—so) — |(io(~'So)l) 


E(-so) 

Uz ^(1 - 4 ^ 5 ) sgn(do(so)) t'(Mo(so)| - E(so)), 


sgn(do(-So))d(ldo(-So)l - E(-so)) 

( 29 ) 
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where magnetic length I = ^^hc/\eB±\ and 65 = B^/B±. We find the following solutions of the gap equation (| 2 ^ 
(solutions in the case of purely magnetic field are considered in Appendix IB II) : 

(i) Ferromagnetic (F) solution: 


Ao = Ao = Aa; = Aj; = 0 , A^ = =Fy(uo + - 2 u_l), = ez ± + Uz + 2uj_). ( 30 ) 

This solution exists for jA^I < |Az|. Using Eqs.dl]) and (| 27 )) . it is easy to check that, for 65 —>• 0 , this solution is 
characterized by the unique order parameter ('I'Ictz'I') which defines uniform magnetization. According to Eq.(l 7 |), 
the corresponding term A^cr^y^y^ in the inverse full propagator describes the Haldane-type mass antisymmetric 
in spin. 

(ii) Antiferromagnetic (AE) solution: 


Ao = Ao = Aj, = Aj, = 0, Az = ±^{uo + Uz-2uj_), A^ = ez T ^{uo + Uz + 2u±). (31) 

This solution exists for jAzj > [Azj. If we neglect Az compared to A^, the unique order parameter characterizing 
this solution is (^ct^'I') which defines staggered magnetization. Obviously, the term A^cr^ in the inverse full 
propagator © corresponds to the Dirac-type mass antisymmetric in spin. 

(hi) Canted antiferromagnetic (CAE) solution: 


Aq = Ao = Aa; = 0, Aj, = (^1 “ cc^^) “ 2 u_l) sin6>, 

A^ = ^^^^{uo + Uz- 2 uj_), Az = - „ iuo + Uz- 2 uj_), cosO = ( 32 ) 

2 2 cos 6^ 2uz 

This solution exists when A^, is real. For 65 —^ 0 , this solution is characterized by two nonzero order parameters. 
They are the staggered magnetization in the x direction (’FtTa,'!') and the uniform magnetization ('I'lcrj'I') in the 
z direction. It is important that the order parameter (tE'cra;^') can not be transformed into (^cr2 4 ') by means of 
a rotation in spin space without inducing the uniform magnetization in the y direction. 

(iv) Charge density wave (CDW) solution: 


Aa; = Aa; = A2 = 0, Az=€z, Aq = ± - (uq - - 2m_l) , Aq = {uq + Uz + 2u±). ( 33 ) 

For 65 —> 0 , the order parameter which characterizes this solution is (thdt). The corresponding term Aq in the 
inverse full propagator 0 describes the Dirac mass. 

It is easy to check that for 65 0 pseudomagnetic field induces additional staggered magnetization (^cr^vh) ^ 65 in 

the F solution, and vice versa, additional uniform magnetization ('I'Vz^') ^ &5 in the AF solution. As to the CAF 
solution, pseudomagnetic field produces additional staggered magnetization in the z direction. Pseudomagnetic field 
leads also to the generation of the Haldane mass in the CDW solution in addition to the Dirac mass. 

The carrier density for strained graphene in the LLL approximation for > jSsI is given by 

P = 2^ ~ sgn(/ri(so)) 0 (|ml(so)I - E{so)), ( 34 ) 

So 


where y-L^SQ) = — (sqAq + Aq) and E{sq) = y (soAa, -I- A^Y + (sqA^ + A^)^. One can check that the carrier density 
p equals zero for the F, AF, and CAF solutions, whereas for the CDW state p is nonzero, thus prohibiting it. It is easy 
to show that for Hs —>• 0 the density pcDW ~ B5. Consequently, for B5 = 0 , the CDW solution becomes admissible 
and coincides with the corresponding solution in the case of purely magnetic field obtained in Appendix IB II Clearly, 
the F, AF, CAF solutions for B^ = 0 also reduce to those found in Appendix IB II 









B. |B_l| < IBsl 


Similarly to the previous subsection it is more convenient in the case jbsl > 1 to use the following parameters rather 
the coupling constants Gint^ gzz, and g±±: 


U05 


Gint _ g zz 

27rZ| ’ 27r/| 


U5 


ff-L-L 

27:ir 


(35) 


where = s/hc/leB^l and for specificity we take eB^ > 0. The gap equation in the LLL approximation is obtained 
from Ea. (E51) by replacing the full fermion Green’s function with and multiplying both sides by the LLL 

projector (1 —so)/2- This time twelve parameters and A^,Ai, enter in combinations (/ij/+ Aj^) + si 2 (/ij/— Aj^). 

Therefore, we have only six independent variables + A^ and fli, — A^. Without loss of generality we can put 
/^j/ = Ml/ = 0 so that mjy = 0, and we are left only with parameters d^. 

By making use of = 0 and constricting the propagator on the Sq = — subspace, we obtain the following equation 
for parameters d^,: 


— du{si2)(Jv — e-z(^z = — 


(1 + Sl2^5 ^)2(^05 + Uzb) 


di{si 2 )<^i 


— (1 — Si 2 b^ ^) 2 m 5 


E{si 2 ) 

di{ — Si2)(Ti 


0{E{ai2) - |do(si2)|) + sgn(do(si2))6'(|do(si2)| - £'(512)) 


£(—S12) 

+ *12^^^) Sgn(do(s'i2)) ^(Mo(s'i2)I - £(s'i2))i 


0(£(—S12) — |do(~si2)|) + sgn((io(—si2))6*(|do(~si2)| — £(—512)) 

(36) 


where d^ = A^ — si 2 Ay and £( 512 ) = + 

We find the following solutions of the above gap equation (solutions in the case of purely pseudomagnetic field are 
considered in Appendix IB 21) ; 

(i) F solution: 


Aq — Aq — Ax — Ax — 0, Az — —ez T + ^25 + dws), Aj, — ±-(uo5 + Uz5 — dus). (37) 

ZO 5 z 


This solution exists for jAzj < lAzj. 
(ii) AF solution: 


Ao = Ao = Aj, = Aa; = 0 , A2 =-ez ± ^(uo5 + U25 + 4^5), A^ = =F 7 r|-(uo 5 + - dus). ( 38 ) 

z Zb^ 

This solution exists for jAzj > lAzj. 

(hi) CAF solution: 


Aq — Aq — Ax — 0, 

Az = ^(no 5 T Uz 5 


Aa, = 

dlts), A 2 






2652 : 


(uo5 + Uz5 


2(^05 + Uz5 
dus), Z = 


- dus), 


£z 

4m5’ 


(39) 


This solution exists only when A^, is real. According to the analysis performed in Appendix IB 21 the CAF 
solution is absent in the case where only pseudomagnetic field is present. Therefore, the question arises what 
happens with the CAF solution found here as £|| =0 and B± —>■ 0. Obviously, A^ vanishes in this limit because 
z is proportional to ez and tends to zero. On the other hand, b^z does not depend on B± and, therefore, is 
not sensitive to the limit B± —>■ 0. Furthermore, the Hamiltonian of the system in the absence of magnetic 
field is invariant with respect to SU{2) spin rotations because the Zeeman term vanishes in this case. Then 
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solutions with different components (i = x, y, z) in A = /S.xUx + A^cTy + A^cr^ in ansatz 0 but with the 
same y^A^ + A^ + A^ are physically equivalent because they can be easily connected by means of appropriate 

SU{2) spin rotations. It is not difficult to check that i/Al + Af for the CAF solution (15^ equals |A 2 | of the 
F solution (l37l) as well as IA 2 I of the F solution (IBllI) in the case where only pseudomagnetic field is present. 
Thus, the F and CAF solutions found in this subsection are physically equivalent if magnetic field is absent. 

(iv) QAH solution: 


Aa, — Ax — Ajs — 0, A 2 — 62 , Aq — (^05 3'Ujs 5 dlis), Aq — i ^ (^05 ’^2,5 4zi5). (^0) 

ZO 5 z 


Using Ea. p7|) . we find that the carrier density for strained graphene in the LLL approximation equals 

1 


P = 


2 TTli 


■ + *12^5 Sgn(yL(si2)) 0 (|ml(s 12)| “ ^(siz)) , 


(41) 


where ^^( 512 ) = Aq — S 12 A 0 and E{si 2 ) = y (A^ — si 2 Ax)^ + (A^ — si 2 A^) 2 . For the F, AF, and CAF solutions, 
the carrier density equals zero while it is nonzero for the QAH solution. Thus, the QAH state is not realized for 
B± 7 ^ 0. For B± —>• 0, 


PQAH — 


sgnAp 

irhc 


eB±. 


(42) 


Therefore, for B± = 0, the QAH solution becomes admissible and coincides with the corresponding solution in 
Appendix IB 21 Other solutions also reduce to those in Appendix IB 21 if i?|| = 0 and B± = 0. 


V. PHASE DIAGRAM OF = 0 QH STATES IN MAGNETIC AND PSEUDOMAGNETIC FIELDS 


We found above several solutions of the gap equation for quasiparticles in graphene in magnetic and pseudomagnetic 
fields. In order to determine which of these solutions is the ground state, we should calculate their energy densities 
and then find out the phase diagram of the system. We used the Baym-Kadanoff formalism^ in order to calculate 
the energy density H of the system: 


n= - 


1 


|eHj_ - soSi 2 e 55 | 


Airhc 

SQS 12 

1 + si2Sgn(eHj_ - soSi2eH5) 


E + ez 


sodz - ruz 
E 


6 {E - |yi,|) + Mo - somol - E) 


(43) 


where mo,z, do,z,/ii, and E are functions of discrete variables so,si 2 (see Eqs. dsn]),® and definitions of mo,z,(io.z 
after Ea. (fT51) l. The derivation of the energy density H is given in Appendix ICl and we retained in Ea. (H51) only the 
contribution due to the lowest Landau level. Therefore, our analysis is valid when the magnitude of dynamically 
generated gaps is much less than the maximum of Landau gaps ^J2hvp\eB±]/c or yj2hv%\eB^\/c. 

By making use of the energy density (H5]) and the solutions found in Sec lIVI we easily calculate the following energy 
densities for the solutions in the case Ms] < 1 : 


Hp = 


M-l| T ez 

ttP 


— C'<, Qaf = 


\u±\bl ± b^ez 
ttP 


— C'<, V,CAF = — 


4:\u±\ttP 


-C<, C< = 


(1 


bl){uo + Uz- 2u±) 




(44) 


where signs in the expressions for H f and Haf correlate with the corresponding ones in solutions (l30)) and (|M]) . These 
energy densities in the nonstrained limit equal up to a constant to the corresponding energies found in Refs. [2^124 1. 
For Msl > 1, we have 


2M565 ^ ± ezb, 


= — 


'5^ 


nq 


— C'>, Haf = — 


2u5 =F ez 


-C>, U 


>, ^^CAF — 


Su^ttP 


-C'>,C> = 


(I + 65 ^)(U05 + Uz5 - 4 ^ 5 ) 


( 45 ) 
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FIG. 1: (Color online) The phase diagram of the i/ = 0 QH states in graphene in the plane (611,65) for g±± > 0 . The blue 
color represents the CAF state, green - AF state, and yellow - F state. This diagram is obtained for B± = 20 T, u± = — 13 K, 

2|e|(;xx/(/rs7r6c) = 2.8, and ez = 13.4y^l + 6^| K. 



V 


i 

\ 


/ CAF 



F 

u\ 


\ 


FIG. 2 : (Color online) The phase diagram of the v — 0 QH states in graphene in the plane (6||, 65) for —^s7r6c/(2|e|) < g±± < 0 
(left panel) and g±± < —/rB'7i'6c/(2|e|) (right panel). The blue color represents the CAF state, green - AF state, and yellow - 

F state. These diagrams similarly to Fig[T|are obtained for B± = 20 T, u± = — 13 K, and ez = 13 . 4 y^l + 6^| K. For the left 
panel, 2|e|gxx/(/rB7i'6c) = —0.8, whereas 2|e|(;xx/(/rs7r6c) = —2.8 for the right panel. 


where signs in the expressions for Q p and correlate with the corresponding ones in solutions (1571) and (1551) . The 
terms C'< and C'> describe a common shift in energy density for three solutions, therefore, they are not important for 
determining the ground state. Hence the energy density depends only on two essential four-fermion coupling constants 
g±z and g±± (related to u± and U5, respectively). 

Using the energy densities Eos. (1551) and (1551) . we find the phase diagram of the iz = 0 QH states in graphene in 
constant magnetic and pseudomagnetic fields. More precisely, we fix the value of perpendicular magnetic held B± 
and study the phase diagram in the plane ( 611 , 65 ), where 6 || = B\\/B± and 65 = B^/Bp. In our analysis, following 

Kharitonov^^, we assume that u± < 0. Then for 2 |mx| > ez (recall that ez = /rs|Hx |+ 6 ^|), we hnd that the 
phase diagram depends on coupling constant g±±. 

The simplest phase diagram is realized for positive g±±, which we plot in FigUJ The CAF solution is found to be 
the ground state of the system in the region dehned by I65I < ez{h\\) /2\u±\ < 1, where this solution exists. The side 
borders of the region with the CAF solution are determined from the condition I 65 I = cosO = ez(6||)/(2|'Ux|). This 
makes the central blue regions, where the CAF state is realized, look like biconcave lens. The top and bottom borders 
of the region with the CAF solution are determined from the condition cosO = ez{b\\)/{2\u±\) = 1. 
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The ferromagnetic solution with the minus sign in Ea. (l44l) is the ground state for |&||| > and I 65 I < 1, where 
is determined by ez{b^^) = 2|ux|, and the AF solution is the most preferable solution for the regions I 65 I > 1 and 
| 6 ||| < b^^, 1 > 1 65 1 > e 2 ( 6 ||)/ 2 |Mx|- As to the AF solution, plus and minus sign for it in Fq.(j44]) is preferable for 
65 < 0 and 65 > 0, respectively. The regions of the CAF, F, and AF states are shown in Fig[T] by blue, yellow, and 
green colors, respectively. 

The phase diagram of the system is somewhat more complicated for negative g±±. There are two qualitatively 
different cases, namely, —fiB'!^fic/(2\e\) < g±± < 0 and g±± < —^B 7 r/ic/( 2 |e|). The corresponding phase diagrams 
are plotted in the left and right panels of Fig 121 respectively, where like in Fig [Tithe blue color represents the CAF 
state, yellow - F state, and green - AF state. Note the symmetry of the phase diagrams with respect to the change 
65 —>• —65 and 6 || —>■ —1>||. The main difference of the phase diagram in the left panel of Fig |2l compared to that in 
FigHlis the appearance of two additional regions, where the CAF state is realized as the ground state of the system. 
The reason for this is that the terms with U 5 in the energy densities of the AF and CAF solutions (l45ll for I 65 I > 1 
are of different sign. Therefore, for g±± > 0, the energy density of the AF state is smaller than that of the CAF 
state. The situation changes for g±± < 0, where the energy density of the CAF state is smaller than that of the AF 
state for sufficiently small 6 ||. As 6 || increases, the CAF state ceases to exist for ez > | 4 it 5 | and transforms into the 
AF state. The side borders of two regions with the CAF solution in the left panel in Figlll are determined by the 
condition (^s7r^c/(2|e||gx_L |))= l&sl- 

For g±± < —/rB 7 r/ic/( 2 |e|), the phase diagram becomes even more complicated. Though in this case for I& 5 I > 1 
the CAF solution is also preferable in every point where it exists, this region is now separated in four detached parts. 

The CAF-AF borders are analogously determined by the condition {2\e\\g±±\))^ 1 + = l^sj. The CAF-F 

borders are constant in 6|| and determined from the condition {2\e\\g^_±_\/ {^bt^^ic)) = . The F solution is 

more preferable than the AF solution for small 6 || and I 65 I > 1 as is seen from the energy densities in Fa. (H51) . 

It is instructive to compare the obtained results with those in Refs. [^l^ where purely magnetic field was considered. 
The case of purely magnetic field corresponds to the line 65 = 0 in our phase diagrams. According to Figsjl] and |2l 
there is the continuous phase transition from the CAF state to the F state as parallel magnetic field B \| increases that 
agrees with Kharitonov’s findings. 


VI. SUMMARY AND DISCUSSIONS 

In the present work, we studied the ^ = 0 quantum Hall states in strained monolayer graphene under tilted external 
magnetic field in the presence of local Coulomb and other local four-fermion interactions, which, in general, break 
the approximate spin-valley SU(4) symmetry of the low-energy electron Hamiltonian of graphene. Solving the gap 
equation in the LLL approximation, we found a rich phase diagram where different QH states, the ferromagnetic, 
antiferromagnetic, and canted antiferromagnetic, compete when we change parallel magnetic field i?|| and pseudo- 
magnetic field i ?5 (we keep perpendicular magnetic field Hx fixed and rather large). For zero strain, these states are 
characterized by the nontrivial order parameters ^cr^'I', and and H'icr^fl', respectively. In the presence 

of pseudomagnetic field these order parameters gain an admixture of uniform staggered magnetization or uniform 
magnetization in the 2 direction. 

Assuming that the strength of the Coulomb interaction is larger than all other local interactions, we found an 
essential dependence of the phase diagram only on two four-fermion couplings g±z and gxx, which appear in the 
low-energy effective Hamiltonian ([2|). Our main results are accumulated in Figs |T] and |21 For negative coupling g±z 
we found that the canted antiferromagnetic state is always preferable in the center of all figures, i.e., for not too large 
values of i?|| and H 5 . When i3|| and increase, the ferromagnetic, antiferromagnetic, and canted antiferromagnetic 
states are realized as the ground state of the system depending on the values of i?||, H 5 . 

The particular case of purely magnetic field corresponds to the line H 5 = 0 in our Figs. I and 2 where there is a phase 
transition from the CAF state to the F state as parallel magnetic field i?|| increases that agrees with Kharitonov's 
findingSi ^^i^^ . In the case of purely pseudomagnetic field considered in Appendix IB 21 we find that the F and QAH 
solutions for g±± < 0 have the lowest and equal energy densities. Since these solutions are described by Haldane-type 
masses, our results agree with the findings in Ref.[20|, where it was shown that pseudomagnetic field catalyses the 
generation of the Haldane mass. On the hand, for g±± > 0, the AF state with Dirac-type mass is the most preferable. 
These results and those obtained in Sec fVl show that the structure of the phase diagrams in Figs. [T]and|2]is sensitive 
to the sign and the strength of four-fermion coupling g±± as well as coupling g±z- 

In our analysis, we ignored boundaries of graphene. However, it is known that they may play an essential role in 
strained graphene. A novel magnetic ground state, where the Neel and ferromagnetic orders coexist, was reported for 
the Hubbard Hamiltonian in strained graphene in Ref.ji^. Whereas the Neel order takes the same sign through the 
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entire system, the magnetization at the boundary takes the opposite sign from the bulk. Since the total magnetization 
vanishes, the magnetic ground state is edge-compensated antiferromagnet. 

In conclusion, strained graphene iir a magnetic field provides a uirique opportunity to observe various symmetry 
breaking phases. Further progress in achieving strain induced pseudomagnetic fields, especially uniform pseudomag- 
netic fields, might allow one to probe experimentally the obtained phase diagram. As for the future, it would be 
interesting to extend the results of the present paper beyond the neutral point with the filling factor 1 ^ = 0 and 
describe the quantum Hall states with other filling factors. 
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Appendix A: Quasiparticle propagator: Expansion over LLs 


For the Fourier transform in time of the full propagator G{u, u') in Eq. (|S]) we can write 

G(a;,r,r') = i(r| [(w-|-771)7° — u_F(7r • 7) — d] ^ |r 0 = *(i’l [(“^ + ~• 7)-I-d] 

X [((w -I- 777)7° ~ • 7) — d) ((w -I- 777)7° ~ vf{t^ ■ 7) + d)] ^ |r') 


= i[W - VF{T^r ■ 7)] (r| {M - VpTv'^ - 7 eHj_(?i 7 ;|/c) 7 ^ 7 ^) ^ |r'). 


(Al) 


where m = fi + 77°7^7^/i, d = A — 77°7^7^ A, the matrices n, p,, A, A are defined after Eq. o, and matrices W and 
M are 


with 


W = (w -I- 777 ) 7 ° -I-d, 777 = mjyCTjy, d = di,au^ f = 0 ,x,z, 
M = c + Ciai, i = x,y,z, 


c = + 2a;777o -I- — d^, Cx = 2 {ujmx + rnorrix — dodx), 

Cy = 2i-j°{dxmz - dzUix), Cz = 2{umiz + rnoiriz - dodz). 


(A2) 

(A3) 


(A4) 


Our aim is to find an expression for the propagator m as an expansion over LLs (we follow below the consideration 
in Appendix A in Ref.1. The operator tt^ has well known eigenvalues {2n + l)h\eBj_\/c with 77 = 0,1, 2,... and 
its normalized wave functions in the Landau gauge A = (0,B±x) are 


i’np(r) = 


\/2ttI y/TBlUjF 


H, 


{j+pl) e"^ 




(AS) 


where Hn{x) are the Hermite polynomials and I = y/hc/leBxl is the magnetic length. These wave functions satisfy 
the conditions of normalizability 


and completeness 


J d^r-ilj*p{r)ipn'p'{r) = Snn'S{p-p'), 

^ 00 

00 p 

/ °^PCp(r)^np(r') =d(r-r'). 


(A6) 


(A7) 


Using the spectral expansion of the unit operator (IA7I1 . we obtain 

(r| {M - - ieBF{hvl/c)'^^'^'^)~^ |r') = exp 


{r-r^ _^{x + x'){y-y') 


E 

n—0 


M. — (277 -I- \){hv‘]r/c)\eB f \ — i{hvp/c)eBj_^^^'^ 


Lr, 


{r-v'f 

2P 


2 P 


(AS) 
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where we integrated over p by using 


OO 

J e-^"H^{x + y)Hr,{x + z)dx = 2^TT^/^m\z^-^L^'^{-2yz), 


(A9) 


assuming m < n. Here L“ are the generalized Laguerre polynomials, and L„ = L®. The matrix i{hvp/c)eB 
has eigenvalues ±?it>|,|eHj_|/c = ±e^/2. Therefore, we have 


LuiO 


P-L^iO , P+LniO 


Ai — {2n + l){hvjp/c)\eB±\ — i{hvp/c)eB±j^^'^ Ai — ne^ Ai — {n+\)e^g' 
where ^ = (r — r')^/(2/^) and the projectors P± are 

P± = i [l ± i 7 ^ 72 sign(eHj_)] . 

By redefining n —>■ n — 1 in the second term in Eq. (IA10|) , Eq. (IA8I) can be rewritten as 

(rllX - vW - Y. 

where L_i = 0 by definition and 


n—0 




(AlO) 


(All) 


(A12) 


(A13) 


is the Schwinger phased which appears due to the noncommutative character of magnetic translations^. Since 


( _ig _ 


y-y 

2 P 


, TTye*"^ = e^'^h ( —idy + 


2/2 


(AM) 


propagator (ED can be presented as a product of the phase factor and a translation invariant part G{u!', r — r'), 

G{uj- r, r') = e*‘^('’’''')G(u;; r - r'), (A15) 

where 


G{u}-, r — r') = i 


W — hvF"f^ ( —idx — 


y-y 

2/2 


— hVFj^ —i9y + 


X — X 

~w 




27r/2 


n—Q 


Ai — nep 


The Fourier transform of the translation invariant part of propagator (jAlfip can be evaluated by performing the 
integration over the angle, 


271- 


dee 


ikr cos 9 


= 27 rJo(fcr), 


(A17) 


where Jo{x) is the Bessel function, and then using the following formula: 


j xe Jo(a: 2 /)c/x = 




a 


n+l 


e-^y Lr. 


Py^ 


2a{P — a) 


valid for ?/ > 0 and Rea > 0. We obtain 




Ai — ne% 

n=0 ^ 


with 


D„(w, k) = 2W [P_L„ {2^1^) - P+T„_i (2fc2/2)] + ihvF{^ ■ i)L\_^ {2 kH^) , = 0, 

describing the nth Landau level contribution. 


(A18) 


(A19) 


(A20) 
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Appendix B: Solutions of gap equation in purely magnetic or pseudomagnetic field 


1. Purely magnetic field 


Let us find solutions of the gap equation (1241) in the LLL approximation in the case where only magnetic field is 
present. In this case, the gap equation takes the form: 


- du{so)cr^ + s^ezCTz = - ^{uo + Uz) 


di{so)c 


Eiso) 


'-9{E{sq) - Mo(so)l) + sgn(do(so))6'(Mo(so)| - E{so)) 


+ u± 


di{ 


9{E{—so) — |c?o(—So)l) + sgn(do(—So))^(Mo(~'So)l ~ E{—so)) 


[ Ei-so) 

Uz sgn{do{so))0{\do{so)\ - E{so)) 


(Bl) 


where dv = Ai, + sqAi, and E = \/+ d^. Projecting on the Pauli matrices, we get 


c^o(so) = 2 ('“0 - u^) sgn(do(so)) 6'(Mo(so)| - ^^(so)) - (mu + u^) sgn(do(-so)) 6'(Mo(-so)| - E{-s^)), 


(B2) 


dx{sQ) = i(uo + Uz)^^difA0(^E{so) - Mo(so)l) - u±^^j—^0(E(-so) - |do(-so)|), 
2 E{so) E{-so) 


(B3) 


dz(so) - soez = ^(uo + Uz)^^7^9(E(so) - Mo(so)l) - u±^^7-^^9(E(-so) - |do(-so)|). (B4) 

2 E(so) E(-so) 

These equations define a system of non-linear equations because £1 is a non-linear function of d^ and d^- Since Gint 
approximates the Coulomb interaction, which is the strongest electron-electron interaction in graphene, we assume in 
what follows that uq ^ ez, Uz, u±. The system of equations (IB2p - (IB4l) has the following solutions which are consistent 
with the charge neutrality condition for both ± signs: 

(i) F solution: 


Aq — Aq — Aa, — Aa; — 0, 


A^ = 0, A^ = ±-(uo + -I- 2m_l). 


(ii) AF solution: 


Ao = Ao = Aa; = Aa; = 0, Aj = ±i(uo -I- - 2u_l), Az = ez, 

which exists for |Aa;| > ez that is satisfied automatically since we assumed uq ^ ez,Uz,u±. 

(iii) CAF solution: 


(B5) 


(B6) 


Ao = Aq = 0, Aa; = A^ = 0, cos 0 = - 


ez 


Ax = ±^{uo + Uz - 2u±) sin0, A^ = ^(uq + Uz - 2u±) cos6. 


(B7) 


The CAF solution exists only for ez < 2|m_l|. 

(iv) CDW solution: 


Aq = Aa; = Aa; = A^ = 0, Az=ez, Ao=±-{uo-3uz-2u±). 


(B8) 
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The ferromagnetic and CAF solutions reproduce solutions obtained in Refs. j^l^ . Note that here like in SecfVlwe 
follow Kharitono v^^i^^ and assume that itj_ < 0. In this case, the CAF solution is the ground state of the system for 
sufficiently small parallel magnetic field. 

Since in the present paper we consider the state with zero filling factor, i/ = 0, we need to control the carrier density 
p which according to Eqs. (HH),® in the LLL approximation is given by the expression 

P = ^ (Mi(so)) 9 (lAiL(so)l - E{so )), (B9) 

So 

where ^l(so) = —(so^o + Aq) and E{sq) = (Aj, + soAa;)^ + (Aj + soA2)2. It is easy to check that for all obtained 
solutions the carrier density p vanishes so that the filling z/ = 0 is realized. The ground state is determined by the 
solution with the lowest free energy density. The case of purely magnetic field corresponds to the line Bs = 0 in Figs.I 
and 2 in Sec|Vl One can see that there is a phase transition from the CAF state to the F state as parallel magnetic 
field Bji increases that agrees with the results obtained by Kharitonov in Refs. [2M^ |. 


2. Purely pseudomagnetic field 


In this subsection, we consider solutions of the gap equation (I26p in the case where only pseudomagnetic field is 
present (i?j_ = R|| =0) retaining only the LLL contribution. Without loss of generality we can put = pi, = 0 so 
that rrii, = 0, and we are left only with parameters di^. Moreover, due to the absence of the Zeeman term the gap 
equation possesses the SU(2) spin symmetry, so that we can take dx = 0. Therefore, we get the following equation 
for the parameters do and dz'. 


- dz{si 2 )(Tz - do{si 2 ) 


2 (^05 T ’^zs) 


sgn[dz{si2))9[E{si2) - \do{si2)\)az + sgn(do(si2))6'(|do(si2)| - E{si2)) 


27 x 5 


Sgn{dz{-S12))0{E{-Si2) - |do(-Si 2 )|)cr 2 + Sgn(do(-Sl2))6'(Mo(-Sl2)| - E{-Si 2 )) 


Uz5 sgn(do(s;2))6'(Mo(s'i2)l - E{s[ 2 )), 
^12 


(BIO) 


where do,z = Aq.z - si 2 Ao,z and E = \dz\. 

Assuming that mqs >> UzSjUs, we find the following solutions: 

(i) F solution: 


(ii) AF solution: 


(hi) QAH solution: 


Ao = Ao = 0, A^ = 0, A^ = ±-(uo5 +Mz5 - 4u5); 


Ao = Ao = 0 , Az = ±-{uo5+Uz5 + Az = 0 -, 


Az — Az — 0, Ag — 0, Aq — ±-(uo5 + Uzo — 4 ^ 5 ). 


(BII) 


(BI 2 ) 


(BIS) 


Note that the CAF solution is no longer present due to the absence of the Zeeman term. For the carrier density we 
obtain the following expression: 

P ~ 2t:P 'Y/ ~ si2Ao^ 9 ^|Ao — S 12 A 0 I — | A^ — si2A2|^ . (B14) 

® S12 


One can check that the condition p = 0 is satisfied for all obtained solutions. The ground state of the system is 
determined as the solution with the lowest free energy density (see SeclV|. By making use of the energy density 
(H51) for B± = = 0, we find that the F and QAH solutions for g±± < 0 have the lowest and equal energy 

densities. According to Eas. (IBll|) and (IB13I) . the E and QAH solutions are described by Haldane-type masses. These 
results agree with the findings in Ref. , where it was shown that pseudomagnetic field catalyses the generation of 
the Haldane mass. On the hand, we find that for g±± > 0 the AF state with Dirac-type mass (IB12|) is the most 
preferable. Once again, like in SeclVl these results show the crucial role of the coupling constant g±± in the selection 
of the ground state of the system in a strong pseudomagnetic field. 
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Appendix C: Free-energy density 


The energy density is defined by the following expression = —TjTV, where TV is a space-time volume. By 
utilizing the Baym-Kadanoff-Jackiw-Tomboulis) formalism, we find the effective action T at its extrema in the mean- 
field approximation (for details see Ref. HI): 


r = —iTr 


\nG-^ + ^{S-^G-l) , 


(Cl) 


where the trace, logarithm, and product S ^G are taken in the functional sense, G = diag(G'+, G_), 1 is the unit 
operator in both matrix and coordinate sense and the expressions for the free and full propagator are given by Eqs. 
and 0. Performing the Fourier transform in time, integrating by parts the logarithm term, and omitting the irrelevant 
surface term, we arrive at the expression (for simplicity in this section we put constants fi = c = 1): 


r = 




-w —+ 2 ^^ iw)G{w) - 1) 


(C2) 


with 


dG ^(w) 
dw 


—i^^5(r — r'). 


The multiplier T came from the functional trace in time, and now Tr contains only spatial integration. By substituting 
the expression for the Green function (IA15I) into ft = —T /TV, one can see that Schwinger phase $ goes away and 
after the Fourier transformation we get for the energy density 


ft = i 


dw 

^ j 

^ dw 
^ 47r 


drk 
f (Pk 


tr 


i^°wG{w, k) -H -[-i{{w - ezcrz) 7 ° - vpik ■ 'y)}G{w-, k) - 1] 


(27r): 


■tr{[(r(; -|- + VF[k ■ 'y)]G{w; k) + i}. 


(C3) 


By making use of the explicit form of the propagator, we calculate the following integrals which contribute to the 
energy density: 


d^k 


{2ny 


VF{k ■ ■j)Gs{w,k) = 


n=0 

i\eB\ 

2 tt 


A4 — ne^ 


E, 

n—O 


A4 — ne% ’ 


(C4) 

(C5) 


where in strained case we use the effective magnetic field B = B± — S 0 S 12 B 5 in the Landau energy es and projectors 

P±- 

By dropping an infinite divergent term independent of the physical parameters and normalizing ft by subtracting 
its value at = 0, we obtain the following expression: 


ft = - 


i\eB\ 

47r 


dw 

—— tr y 
27r 

n—0 


{w + ez<Jz)j°W{P^ + P+6{n - 1/2)) -|- ne| 

A4 — ne^ 








w 


i 2 _ 


(C6) 


After integrating this expression over frequency and taking trace, we finally arrive 


at Ea. (H51) . 
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